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We report a search for B decays to selected final states with the meson: —>■ K'^rjc'K'^'K~, 

—>■ K^rj^ui, B^ —>■ K^rjcTi and B^ —^ K^rjcn^. The analysis is based on 772 x 10® BB pairs 
collected at the T(45') resonance with the Belle detector at the KEKB asymmetric-energy e'''e~ 
collider. We set 90% confidence level npper limits on the branching fractions of the studied B 
decay modes, independent of intermediate resonances, in the range (0.6 — 5.3) x 10~^. We also 
search for molecular-state candidates in the D°D*° — D^D^ + and D*°D*° -|- 

combinations, neutral partners of the ^(3900)^ and Z(4020)^, and a poorly understood state 
X(3915) as possible intermediate states in the decay chain, and set 90% confidence level upper limits 
on the product of branching fractions to the mentioned intermediate states and decay branching 
fractions of these states in the range (0.6 — 6.9) x 10“®. 

PACS numbers: 13.25.Gv, 13.25.Hw, 14.40.Pq, 14.40.Rt 


INTRODUCTION 


Many exotic charmonium-like states are observed in 
the mass region above the DD threshold. Decays of 
B mesons provide a fruitful opportunity to study these 
states and to find new ones. For example, the state 
X(3872) was first observed by Belle in exclusive B+ —>■ 
BT+tt+tt” J/i/) decays [i| that was later confirmed by 
CDF 0, D0 0 and BaBar [^. It was also observed in 
the LHCb experiment [^@] in pp collisions and B decays. 
The X(3872) mass is close to the rrijjo -\- TOp,.o thresh¬ 


old, which engendered a hypothesis that this state may 
be a D^D*^ molecule 0. The observation of the decay 
X(3872) —>• 'yJ/ip by BaBar @ and Belle Q established 
the charged parity of X(3872) to be positive. Angular 
analysis of the A(3872) —>■ J/Tr~ decay by LHCb Q 
determined all its quantum numbers: 

If A(3872) is indeed a D^D*^ molecule, there can ex¬ 
ist other “A(3872)-like” molecular states with different 
quantum numbers. Some may reveal themselves in the 
decays to final states containing the rjc meson. For ex¬ 
ample, a D^D*° — D^D*^ combination (denoted here- 
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inafter by Xi(3872)) with quantum numbers 
would have a mass around 3.872 GeV/c^ and would de¬ 
cay to rjcP and rjcOJ. Combinations of + D^D^, 

denoted by X(3730), and D*°D*^ + D*°D*°, denoted 
by X(4014), with quantum numbers = 0++ would 
decay to pcP and ??c7r°. The mass of the X(3730) state 
would be around 2mjjQ = 3.730 GeV/c^ while that of the 
X(4014) state would be near 2m£).o = 4.014 GeV/c^. 

Recently, a new charged state .Z(3900)* was found in 
Y (4260) decays by Belle [l^ and BESIII [nj. Since this 
particle is observed in the decay to it should 

contain at least four quarks. BESIII [l^ reported sub¬ 
sequently an observation of another decay channel of 
the seemingly same state Z(3885)^ —>■ (DD*)^. The 
Z(3900)* was confirmed in the decay to TT^J/ij) by an 
analysis of CLEO-c data [l^ that also reported evidence 
for its neutral isotopic partner Z(3900)°. Another ex¬ 
otic charged state Z(4020)^ was observed by BESIII in 
decays to ir^hc 0 and {D*D*)^ [T^. There are some 
indications from these analyses that the spin and parity 
of the charged states might be = I"*". 

The near-threshold enhancement in the ujJ/ip invariant 
mass distribution named y(3940) was first observed by 
Belle in exclusive B —?> A'wJ/'i/; decays 0. Later, in the 
same decay mode, BaBar discovered Ar(3915) 0L which 
was confirmed by Belle in two-photon production ISj and 
by other BaBar measurements [1^ l20| . The parameters 
of y(3940) are consistent with those of X(39I5), so they 
are considered to be the same particle. The quantum 
numbers of Al(3915) are claimed to be = 0++, but 
its nature is still undetermined, and there are several 
interpretations describing this state 21-^. 

To search for the particles described above, we recon¬ 
struct r]c mesons via the KgK^n^ mode and study the 
following four decays of charged B mesons: 

1. the (tt+tt”) decay mode: B^ —>• K^X —>• 
K^{r]cTT~^TT~), where we look for Xi(3872), 
y(3900)° and ^(4020)°; 

2. the (w) decay mode: B^ — K^X — K^{r]cOj), 
where we look for Xi(3872); 

3. the (ry) decay mode: B^ —>• K^X — K^{r]cr]), 
where we look for X(3730), X(40I4) and X(39I5); 

4. the (7r°) decay mode: B^ — K^X — K^{r]cn°), 
where we look for X(3730), X(40I4) and X(39I5). 


EVENT SELECTION 


The analysis is based on a data sample that contains 
772 X 10® BB pairs, collected with the Belle detector 
at the KEKB asymmetric-energy e+e“ collider 2^, |27| 
operating at the T(45') resonance. 

The Belle detector [2^ (also see detector section 
in [ 2 ^) is a large-solid-angle magnetic spectrometer that 


consists of a silicon vertex detector, a 50-layer central 
drift chamber (CDG) for charged particle tracking and 
specific ionization measurement {dE/dx), an array of 
aerogel threshold Cherenkov counters (ACC), time-of- 
flight scintillation counters (TOE), and an array of 8736 
CsI(Tl) crystals for electromagnetic calorimetry located 
inside a superconducting solenoid coil that provides a 
1.5 T magnetic field. An iron flux return yoke located 
outside the coil is instrumented to detect AT® mesons and 
identify muons. We use a GEANT3-based Monte Carlo 
(MC) simulation to model the response of the detector 
and determine its acceptance [ 3 ^ . 

Charged tracks are selected with requirements based 
on the goodness of fit of the tracks and their impact pa¬ 
rameters relative to the interaction point (IP). We require 
that the polar angle of each track be in the angular range 
(18° —152°) and that the track momentum perpendicular 
to the positron beamline be greater than 100 MeV/c. 

Pions and kaons are distinguished by combining the 
responses of the ACC and the TOE with dE/dx mea¬ 
surements in the CDC to form a likelihood ratio Cijj = 
Cil{Ci + Cj). Here, £i is the likelihood that the parti¬ 
cle is of type i. Kg candidates are reconstructed via the 
7r“''7r“ final state. The tt+tt” invariant mass must lie in 
the range 0.486 GeV/c^ < M(7r+7r“) < 0.510 GeV/c^. 
The flight length of the Kg is required to lie within the 
interval (0.1 — 20) cm. The angle ip between the pion pair 
momentum and the line joining the 7r“'"7r“ the vertex to 
the IP must satisfy cos:p > 0.95. 

The invariant masses of intermediate resonances must 
lie within ranges obtained from signal MC and contain¬ 
ing more than 95% of the signal yield: 2.9254 GeV/c^ 
< M{KgK^TT^) < 3.0454 GeV/c^ for the rjc meson, 
0.758 GeV/c^ < M{Tr+Tr-Tr°) < 0.808 GeV/c^ for the 
oj meson, 0.125 GeV/c^ < < 0.145 GeV/c^ for 

the 7r° meson, 0.528 GeV/c^ < XI{jj) < 0.568 GeV/c^ 
and 0.538 GeV/c^ < M(7r+7r“7r®) < 0.558 GeV/c^ for 
the ry meson. 

B candidates are identified by their center-of-mass 
(c.m.) energy difference AE = (^ ■ Ei) — E^, and the 
beam-constrained mass M^c = sjE'^jc} — Ic, 

where Eb = \/s/2 is the beam energy in the T(4S') c.m. 
frame, and pi and Ei are the c.m. three-momenta and en¬ 
ergies, respectively, of the B candidate decay products. 
We define the signal region by |AE| < 0.02 GeV and 
5.273 GeV/c2 < Mbc < 5.285 GeV/c^; the |AA| range 
extends to 0.04 GeV for some decay channels. 

Although the continuum background (e“'"e“ —^ qq^ 
where q = u,d,s,c) is not dominant, we suppress it us¬ 
ing topological criteria. Since the produced B mesons 
are nearly at rest in the c.m. frame, the signal tends to 
be isotropic, while the continuum qq background tends 
to have a two-jet structure. We use the angle ©thrust be¬ 
tween the thrust axis[? ] of the B candidate and that of 
the rest of the event to discriminate between these two 
cases. The distribution of | cos ©thrust | is strongly peaked 
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near 1 for qq events but nearly uniform for T(4S') —>■ BB 
events; we require | cos ©thrust | < 0 . 8 . In addition, we 
require that the c.m. polar angle Og oi the reconstructed 
B meson satisfies | cos 0^1 < 0 . 8 , where this angle is mea¬ 
sured relative to the z axis that is collinear with the 
positron beam. 

The mean number of multiple B candidates per event 
varies from 1.2 to 2 . 0 , depending on the decay channel. If 
an event has multiple B candidates, we select the candi¬ 
date with the minimum value of the following expressions 
in the order described below: 

1 . Im^fo - M(7r+7r-)|, 

2 . \mr,,- 

3 . \mrf — M(77)| (for ry —>■ 77 in the (7) mode) or 
|m^o — M(77)| (for the (w) mode, for 77 —7r+7r“7r° 
in the (77) mode and for the (tt^) mode). 



FIG. 1. The A.E distribution for the decay —>■ 

In the (tt+tt”) and (tt^) modes, we observe some sig¬ 
nificant signal and so perform a two-dimensional fit of 
the KsKtt invariant mass (a;) and AE (y) distributions: 

f{x,y) = NsU{x)v{y)+Nnon-resViy)+Co + CiX + C 2 y, (2) 


4. maximum difference between z coordinates at the 
closest-distance point to the z axis among each pair 
of charged particles in the signal-i? final state. 

According to MC, this procedure selects the correct B 
candidate with 95% probability. 


RECONSTRUCTION OF FINAL STATES 

R* —>• + hadrons 

We search for molecular-state candidates 

Z(3900)°, Z(4020)°, and A(3915) in the following B me¬ 
son decays: B^ —>■ K^ric'x'^'x~ ^ B^ —>■ K^rjcCo, B^ —>■ 
K^rjcTj^ and B^ —s- K^rjc'x^. To determine the branch¬ 
ing fractions, we perform a binned maximum-likelihood 
fit of the AE distribution that is modelled by a peak¬ 
ing signal and featureless background. For the (tt+tt”), 
(w) and (tt^) decay modes, the signal function is the sum 
of two Gaussians (G) and the background function is a 
linear polynomial; the fitting function is 

f{x) = Ns [Q:G(a;i,CTi) -f (1 - a)G{x 2 ,cr 2 )\ + cq -h cix. 

( 1 ) 

For a given mode, we generate signal MC and fix the 
mean values xi and X 2 , the standard deviations (Ti and 
tT 2 , and the fraction of the first Gaussian a; we also ob¬ 
tain the detection efficiency for the mode. Here and in 
the following the detector resolution is taken into ac¬ 
count. To account for the difference in resolution between 
MC and data, we replace the resolution of the dominant 
second Gaussian in the fit with (T 2 = \/<x 2 + where 
the resolution degradation 5 = (7.1 ± 2.3) MeV/c^ is 
taken from the analysis of the decay B^ K^rjci^S) —>■ 
K^{K^K^'k^) [31|. The AE distribution for the (w) 
mode is shown in figure [TJ 


u(a;) = 6(M^^,r^J 0 G(0,(Tres), (3) 


r aG(yi, CTi) -b (1 - a)G{y2,(T2) 
v{y) = I (4) 

I Glg( 2 /o,ct,P) 


for the ( 7 r''' 7 r ) and (tt®) modes, respectively. The loga¬ 
rithmic Gaussian function is defined as 


G'LGi.xo.a, P) 


P 



{ 1 -P{x-xq)/ a) 


la 


T 

0 


2 


where (Tq = ^{PV‘2\n2) and P is the asym¬ 

metry parameter. The function u(x) characterizes the 
Tjc resonance and is described by the convolution of a 
Breit-Wigner ( 6 ) function and a Gaussian detector res¬ 
olution function with Cres = 6.2 MeV/c^ obtained from 
Ref. jsij. The parameter Anon-res represents the number 
of events that do not contain an intermediate r]c meson, 
but have the same final state. According to the fit, most 
of the events in the AE peak are of that origin. The 
M(KsKtt) and AE distributions for (tt+tt”) and (tt^) 
modes are shown in figures [2] and |3l respectively. 

To improve the AE resolution in the ( 7 ) mode, we 
modify the energy of the 7 candidate decaying into 
photons: E'^ = + p^, where = 547.853 

MeV/c^ 33| is the mass and Prj is the reconstructed mo¬ 
mentum. Since the rj candidate is reconstructed in two 
decay modes, we perform a combined fit of the AE dis¬ 
tribution corresponding to 77 —?> 77 and 77 —>■ 7 r+ 7 r“ 7 r°, 
using the following function: 


fi{x) = NeffSiBi [aiG(xi,i,CTi,j) -b (1 - ai)G(x2^^,a2,^)] + 
co.i + ci,iX, (5) 

where i refers to either rj —> 77 or 77 —^ de¬ 

cay. In particular, H 27 = B(r] —>■ 77 ) and = B(r] —>■ 
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Decay mode 

Fit. function 

Efficiency, % 

Yield 


0 

0.53 + 0.01 

-41 + 27 

B^ —>■ Ar*7c7r’’'7r' 

0 + 0 + 0 

2.84 ± 0.02 

155 + 72 



3.69 + 0.01 

-1.9 + 12.1 

B^ ->■ K^rtcV, 

7 —>■ 77 

7 —>■ 7 r’’' 7 r~ 7 r'' 

© 



3.05 + 0.01 

-14 + 26 

0.69 + 0.01 

-1.8+ 3.4 


TABLE I. Fit results for B decays independent of intermedi¬ 
ate resonances. 


FIG. 2. Projections of the two-dimensional fit in KsKn 
invariant mass (left) and AE (right) for the decay ^ 

K^r]cTr^'n'~. Each projection is plotted for events in the whole 
fitting range of the other projection. 




FIG. 3. Projections of the two-dimensional fit in KsKtt 
invariant mass (left) and AE (right) for the decay B^ 

Each projection is plotted for events in the whole 
fitting range of the other projection. 


7 r+ 7 r“ 7 r°) x B{Tr'^ —^ 77 ). Here, the parameter Ngs is 
the effective number of signal events. To obtain the total 
yield for each decay channel, it should be multiplied by 
the corresponding efficiency £ 2^/377 and rj decay branch¬ 
ing fraction Bi. The combined fit projections of the AE 
distributions for the {rj) mode are shown in figured 




FIG. 4. The combined fit projections of the AE distribntions 
in case of the 7 —>■ 77 (left) and 7 —>■ 7 r"'" 7 r“ 7 r*’ (right) modes 
for the decay B^ —>■ K^rjarj. 


The fit results are summarized in table ID 


Xi(3872), X(3730) and X(4014) 

For the (tt+tt") mode, we perform a maximum- 
likelihood fit of the r]cTT'^Tr~ invariant mass distribution 
that is, again, modeled by a peaking signal and a smooth 
background; the fit function is 


f{x) = Ns [aG{M, CTi) + (1 — a)G(M, 0 - 2 )] + cq + cix. 

( 6 ) 

For each intermediate resonance, we generate signal MC 
that incorporates the corresponding “X(3872)-like” par¬ 
ticle quantum numbers. From the signal-MC fit, we fix 
the mean value M, the standard deviations tri and (T 2 , 
and the Gaussian fraction a; in addition, we obtain the 
detection efficiency for the mode. From the fit, we ob¬ 
tain the signal yield Ng, which is shown in table HIl The 
invariant mass distribution is shown in figure [5] 

(left). 

We validate our procedure by applying it to the de¬ 
cay ^ K^%l}{2S), 'il’(2S) —>■ J/'i/jTT+TT”. This decay is 
similar to the (tt+tt”) decay except that we reconstruct 
the 'ip{2S) meson in place of the Xi(3872) and the J/ij} 
in place of the 7 c- The J/ip meson, like the 7 c, is re¬ 
constructed via the final state. The selection 

criteria for this decay are the same as for the (tt+tt”) 
decay except for the invariant mass of the 
combination; 3.077 GeV/c^ < M{KgK^'!T^) < 3.117 
GeV/c^. The mean number of B candidates per event 
is 1.7. In case of multiple B candidates, the one with 
the minimum differences for \mj^o — M{'k~^'k~)\ and 
\mj/Tp — M{K^K^'k^) \ and the best vertex coordinate is 
chosen. We fit the invariant mass distribution 

and obtain the number of signal events Ng = 20.2 ± 6.5, 
which corresponds to a significance of 3.5 standard devi¬ 
ations (a). The significance is estimated using the value 

of —27/in , where £max {Bq) denotes the likelihood 


value when the yield is allowed to vary (is set to zero). 
The expected number of events estimated u sing the world 
averages of the known branching fractions 33| is 22 ± 4, 
which is consistent with Ng. 

In the analysis of the (w) mode, we use the sum of two 
Gaussians to describe the signal and a threshold square- 
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root function to describe the background: 

f(x) = Ns [aG{M, CTi) + (1 - a)G{M, 02 )] + - ci. 

(7) 

The rjcUJ invariant mass distribution is shown in figure [5] 
(right). 




FIG. 5. The distributions of ricTT'^Tv~ (left) and rjcui (right) 
invariant masses in the search for the Xi(3872). The thresh¬ 
old of the i]cto invariant mass is slightly shifted relative to the 
sum of the uj and r/c nominal masses due to the mass windows 
of these two mesons, as confirmed by MC studies. 

In the X(3730) mass region of the (r?) mode, the fitting 
function is 

fi{x) = NeBEiBi [aiG{Mi, ai^i) + (1 - ai)G{Mi, a 2 ,i)] + 
co,i + ci^iX, (8) 

where A^eff is the effective number of signal events and i 
refers to either ij —>■ 77 or 7 —>■ 7 r+ 7 r“ 7 r° decay. The rjcT] 
invariant mass distributions in the Ai(3730) mass region 
are shown in figure [5] (top). 

In the Ai(4014) mass region of the (ry) mode, the fitting 
function is 


NeS£2jB2jGBG{M2j, (72^, P) + Cq,! + ci^a;, 


f{x) 


^eBESTrB^TT CTr 37^)“!“ 

(1 - q:)G(M3,7, ct2,37i-)] + C 0.2 + ci,2a;. 


(9) 

The rjcT] invariant mass distributions in the X (4014) mass 
region are shown in figure [ 6 ] (bottom). 

For the ( 7 r°) mode, we use 


f{x)=NsGBG{M,a,P)+co + cix. (10) 

The rjcTT^ invariant mass distributions are shown in fig¬ 
ure [T] 

The fit results are summarized in table m 


^(3900)“ and Z(4020)° 

We perform a sequence of binned maximum likelihood 
fits of the r]c'7'^'7~ invariant mass using the convolution 




M(r|er|), GeV/c^ 




MCri^ri), GeVk^ 


FIG. 6 . The combined fit projections of the 7 c 7 invariant mass 
distributions in case of the 7 —>■ 77 (left) and 7 —>■ 

(right) modes corresponding to the search for the X(3730) 
(top) and X(4014) (bottom) resonances. 




FIG. 7. The ticTv^ invariant mass distributions correspond¬ 
ing to the search for the X(3730) (left) and X(4014) (right) 
resonances. 


of a Breit-Wigner and a Gaussian for the signal and a 
linear polynomial for the background: 

f{x) = b{M, T) 0 G(0, CTres) + Cq + CiX. (11) 


The Gaussian models the detector resolution, which is 
assumed to be similar to that obtained in Ref. 31| and 
equal to 9.8 MeV/c^. The Breit-Wigner mass is confined 
to a 20 MeV/c^ window (the so-called mass bin) that is 
scanned in 20 MeV/c^ steps across the range (3.79 — 4.01) 
GeV/c^ for the Z(3900)° and (3.93—4.07) GeV/c^ for the 
Z(4020)°. The width is fixed to the weighted mean of the 
previously measured values (35 MeV/c^ for the Z(3900)° 
and 12 MeV/c^ for the 0(4020)°). The background is de¬ 
scribed by a linear polynomial. The detection efficiency is 
obtained from signal MG: (9.64±0.03)% for the 0(3900)° 
and (10.42 ± 0.03)% for the 0(4020)°. The obtained sig- 
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Decay mode 

Fit. function 

Efficiency, % 

Yield 

Xi(3872) 7c7r+7r' 


7.95 ± 0.02 

17.9 ± 16.5 

Ai(3872) — )• rjcOJ 

0 

1.92 ±0.02 

6.0 ± 12.5 

A(3730) ^ 

7 — >■ 77 

7 — >■ 7 r'*' 7 r~ 7 r° 

m 



6.57 ±0.02 

13.8 ±9.9 

1.18 ± 0.01 

1.4 ± 1.0 

X(3730) 7c7r“ 

uni) 

6.52 ±0.02 

-25.6 ± 10.4 

X(4014) 7 c 7 , 

7 77 

, + - 0 

77 ^ TT TT TT 

® 



7.09 ± 0.02 

8.9 ± 11.0 

1.78 ±0.01 

1.3 ± 1.6 

X(4014) ->■ 7c7r“ 

uni) 

7.55 ±0.02 

-8.1 ± 13.2 


TABLE II. Fit results for the Xi(3872), X(3730) and X(4014) 
resonances. 


nal yield is shown in figure [S] The examples of the fit 
within the mass bin containing 3.9 GeV/c^ for Z(3900)° 
and 4.02 GeV/c^ for Z(4020)° are shown in figurejS! 




FIG. 8 . Dependence of the signal yield of the Z(3900)° (left) 
and Z(4020)° (right) on the mass bin. The mass bin is a 20 
MeV/c^ window to which the mass is confined and scanned 
in 20 MeV/c^ steps across the fit range. 




FIG. 9. Examples of the fit within the mass bin containing 
3.9 GeV/c^ for Z(3900)° (left) and 4.02 GeV/c^ for Z{4020)° 
(right). 


X(3915) 

For the rjcT] invariant mass distribution, we perform a 
combined fit of two decay modes of the rj meson: 

fi{x) = [b{Mi,ri) 0 G(0,a,)] + 

Co,2 “1“ 

where is the effective number of signal events and i 
refers to either 77 — 77 or 77 — )> 7 r+ 7 r“ 7 r° decay. The detec¬ 
tor resolution ai is obtained from signal MG, taking into 
account the resolution degradation and is equal to 13.6 
MeV/c^ for 77 —77 and 12.3 MeV/c^ for 77 —>■ 7 r+ 7 r“ 7 r°. 
The corresponding 772/7 invariant mass distribution is 
shown in figure [TUI 




FIG. 10. The combined fit projections of the rjcrj invariant 
mass distributions in case of the 77 —>■ 77 (left) and 7 —>■ 
7 r"*’ 7 r“ 7 r'^ (right) modes corresponding to the search for the 
X(3915) resonance. 

We fit the 77 c 7 r° invariant mass distribution with the 
function in Eq. (HU. The detector resolution is obtained 
from signal MG, taking into account the resolution degra¬ 
dation, and is equal to 15.7 MeV/c^. The corresponding 
TjcTT^ invariant mass distribution is shown in figure 1111 



FIG. 11. The 7 c77° invariant mass distribution corresponding 
to the search for the A(3915) resonance. 

The fit results are summarized in table 11111 


SYSTEMATIC UNCERTAINTIES 


Systematic uncertainties are categorized as follows: 





























































































































































Decay mode 

Fit. function 

Efficiency, % 

Yield 

A(3915) —>• 77c 77, 

77 — >■ 77 
rj —>■ 

di) 



6.60 ± 0.02 

-7.4 ± 14.5 

1.64 ±0.01 

- 1 . 1 ± 2.1 

X(3915) - 7 - 77 c 7 r“ 

m 

6.88 ± 0.02 

-4.3 ± 18.1 


TABLE III. Fit results for the A(3915) resonance. 

1. Additive systematic uncertainties affect the num¬ 
ber of signal events and are estimated by the varia¬ 
tion of the fit conditions. They are displayed as the 
numbers of events in tables Hvl and Ivl and arise from 
the sources listed below. (i)To obtain the error re¬ 
lated to the resolution degradation, we vary the cor¬ 
rected variance of the fitting function within its sta¬ 
tistical uncertainty, (ii) We assume that the com¬ 
binatorial background can be parameterized with a 
first-order polynomial. To obtain the background 
shape uncertainty, we describe the background by a 
second-order polynomial and compare the results. 
For the (w) mode, we change the square root to the 
fourth root, (in) To estimate the systematic error 
associated with the selection criteria, we relax the 
criteria on Aif, Mbc and the invariant masses of 77 c, 
Lo and 77 by 50%. (iv) We vary the bin size between 
2.5 to 7.5 MeV/c^ and determine the corresponding 
systematic error. 

2. Multiplicative systematic uncertainties affect the 
product branching fractions. They are displayed in 
percent in table fVTl and arise from the sources listed 
below, (i) The number of BB pairs is calculated 
from the difference between the number of hadronic 
events on resonance and the scaled number of those 
off resonance. The systematic error is dominated 
by the uncertainty in the scale factor and is equal to 
^1.4% [^. (ii) The uncertainties on the uj, 7r°, 77, 
77c and Kg decay branching fractions are taken from 
Ref. [ 3 ^. (Hi) The statistical error of the efficiency 
determined by the signal MC is also taken as a sys¬ 
tematic uncertainty. For Z(3900)° and Z(4020)° 
decays, we take into account the efficiency varia¬ 
tion in the ricP and 770/0 decay modes. For the de¬ 
cays B^ —>■ Rr^77c7r+7r“ and B^ —?> iF±77 c7r°, we 
assume alternative decay modes containing inter¬ 
mediate iF*°, AT*(1410)^ and p mesons, and take 
into account the difference of the MC detection effi¬ 
ciency as the corresponding systematic uncertainty. 

(iv) An analysis of the charged track reconstruc¬ 
tion uncertainty as a function of particle momen¬ 
tum gives an estimate of 0.34% per charged track. 

(v) To determine the errors due to K and tt meson 

identification, data from the analysis of the process 
£)*+ —7> iD°7r+ followed by the decay —>■ K~Tr'^ 

are used. The uncertainty in identification is 
0.8% per K meson and the corresponding value for 


TT^ identification is 0.5% per tt meson, (vi) Esti¬ 
mation of the contribution of the 77 and 7r° recon¬ 
struction uncertainty is carried out using the com¬ 
parison of the number of reconstructed 77 —>■ 37r*^ 
and ?7 —>■ 77 events. Such an estimate gives 2% 
per 77 and tt^ meson, (vii) The contribution of the 
Kg reconstruction uncertainty is estimated to be 
4.4% [^. (via) We also take into account the devi¬ 
ation of MC from the data by applying a correction 
to the efficiency: EData/eMC is 0.9996 for each kaon 
and 0.9756 for each pion. 

For the (77) mode, the rj meson is reconstructed in two 
decay modes: 77 — >■ 77 and 77 — >■ tt^tt~tt^ . To estimate the 
systematic uncertainty, we compare the systematic errors 
given by each 77 decay mode and take the maximum as 
the uncertainty for the reconstruction. 

For the Z(3900)° and Z(4020)° decays to 77 c 7 r+ 7 r“, 
the additive systematic uncertainty is assumed to be the 
same as in Ai(3872) —>■ 77 c 7 r+ 7 r“. In addition, we vary 
the resonance width in the (15 — 65) MeV/c^ interval 
for the Z(3900)° and (2.5 — 27.5) MeV/c^ interval for the 
Z’(4020)°. The intervals are chosen according to the vari¬ 
ation of the previously measured Z(3900)^ and Z(4020)^ 
width values. The total additive systematic uncertainty 
is 28.5 events for the Z(3900)° and 25.9 events for the 
Z(4020)°. 

For the A(3915) decays, the systematic uncertainty 
is taken from the similar A(4014) decays (see tables IVl 
and IVII) . 


Source 

TjcTT^TT 

7?cW 


77c7r" 

Background parameterization 

1 

44 

2687 

2 

Selection criteria 

< 1 

< 1 

1695 

33 

Bin size 

18 

2 

430 

9 

Total (events) 

18 

44 

3206 

34 


TABLE IV. Additive systematic uncertainties (in events) for 
B decays without intermediate resonances. For the pcP mode, 
the uncertainty corresponds to the effective number of events 

Aeff. 


Source 

rjcTT^TT 

77 c 07 

VcV 

Vc 


Mass, MeV/c'‘ 

3872 

3872 

3730 

4014 

3730 

4014 

Resolution degradation 

1.2 

< 1 

68 

28 

1.4 

0.2 

Background parameterization 

5.8 

3.6 

18 

8 

0.8 

0.3 

Selection criteria 

23.9 

5.4 

293 

280 

5.2 

9.3 

Bin size 

1.2 

7.7 

30 

71 

2.4 

4.4 

Total (events) 

24.7 

10.1 

303 

290 

5.9 

10.3 


TABLE V. Additive systematic uncertainties (in events) for 
B decays containing the A(3872)-like particles. For the pcP 
mode, the uncertainty corresponds to the effective number of 
events Nes. 
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Source 

77c7r’’'7r 

ricOJ 


77c7r" 

Number of BB pairs 

1.4 

1.4 

1.4 

1.4 

B{u} —>■ 7r+7r“7r°) 

— 

0.8 

— 

— 

B{tt° — >• 77 ) 

— 

< 0.1 

< 0.1 

< 0.1 

B{ri 77 ) 

— 

— 

0.5 

— 

B(ri —>■ 7r'''7r“7r°) 

— 

— 

1.2 

— 

Bin, ^ 

6.8 

6.8 

6.8 

6.8 

B{K% TT+TT") 

0.1 

0.1 

0.1 

0.1 

MC detection efficiency 





no resonance 

35.8 

2.4 

1.3 

19.5 

Js:(3872)-like 

0.3 

0.5 

0.7 

0.3 

Z(3900)°/Z(4020)° 

13.3/4.4 

— 

— 

— 

Track reconstruction 

1.7 

1.7 

1.7 

1.0 

identification 

1.6 

1.6 

1.6 

1.6 

TT^ identification 

1.5 

1.5 

1.5 

0.5 

77 reconstruction 

— 

— 

2.0 

— 

7 r° reconstruction 

— 

2.0 

2.0 

2.0 

K% reconstruction 

4.4 

4.4 

4.4 

4.4 

Total (%) 

no resonance 

Js:(3872)-like 

Z(3900)“/Z(4020)“ 

36.8 

8.7 

15.9/9.7 

9.3 

9.0 

9.3 

9.2 

21.3 

8.7 


TABLE VI. Multiplicative systematic uncertainties (in %). 


RESULTS 


Decay mode 

Upper limit (90% C.L.) 

— >■ K^rjcTr'^Tf~ 

B± ^ K^ycy, 

3.9 X lO""* 

5.3 X 10-'‘ 

2.2 X lO-'^ 

6.2 X 10“® 


TABLE VII. Results of branching fraction measurements for 
the B decays without an intermediate resonance. 


Resonance 

Decay mode 

Upper limit (90% C.L.) 

Xi(3872) 

ricTT'^TT~ 

3.0 X 10"® 


VciO 

6.9 X 10"® 

X(3730) 

VcV 

4.6 X 10"® 


r]cTv° 

5.7 X 10"® 

X(4014) 

VcV 

3.9 X 10"® 


0 

r]c'K 

1.2 X 10"® 

%(3900)" 


4.7 X 10"® 

%(4020)'^ 


1.6 X 10"® 

^■(3915) 

VcV 

3.3 X 10"® 


0 

r]c'K 

1.8 X 10"® 


TABLE VIII. Results of branching fraction measurements for 
the B decays containing an intermediate exotic resonance. 
For .^(3900)'^ and Z{4020)^ resonances the results are shown 
under the assumption that the masses are close to those of 
their charged partners. 


Our data sample does not permit us to measure the 
branching products of production and decay of the states 
listed above nor the B decay branching fractions so we 
set upper limits instead, taking into account the statis¬ 
tical and systematic uncertainties. The expression for 
the required branching value is z = Ns/^NssBi), where 
Ng is the yield, Nb is the number of BB pairs, e is the 
detection efficiency, and Bi is the product of the branch¬ 
ing fractions of the intermediate resonances. This ex¬ 
pression can be written as z = xy, where x = Ng and 
V — The distribution of the numerator Ng is 

assumed to be Gaussian with mean y, and standard de¬ 
viation Ox = -\Jcr^tat + f^add syst' The distribution of the 
inverse of the denominator is also assumed to be Gaus¬ 
sian with mean v = \/{NBsBi) and standard deviation 
CTy = (Tmuit.syst*^- Thus, the distribution of z can be writ¬ 
ten in the following way: 

F{z)= / G{x,y,cFx)G{-,v,ayv) — dx. (13) 

7-no X \x\ 


The 90% confidence level (C.L.) upper limit C/ on z is 
defined by 


lo Fi.z)dz 
I^F{z)dz 


= 0.9. 


(14) 


Upper limits on the branching fractions and products 
for all the studied decay modes are shown in tables IVIII 
and I VIIH 


Figure [12] shows the dependence on the mass bin of the 
upper limit for the branching product of the Z(3900)° 
and Z(4020)° production and decay: no significant sig¬ 
nal is seen in any of the invariant mass bins. For the 
Z(3900)° resonance, we set upper limits in the range 
(1.8—4.7) X10“^ for the mass region (3.79—4.01) GeV/c^. 
For the Z(4020)° resonance, we set upper limits in the 
range (1.6 — 3.7) x 10“® for the mass region (3.93 — 4.07) 
GeV/c^. If we assume that the Z(3900)° and Z(4020)° 
masses are close to those of their charged partners, we ob¬ 
tain the upper limits on the product branching fractions 
shown in table IVIIII 




FIG. 12. Dependence on the mass bin of the upper limit for 
the Z(3900)'^ (left) and Z(4020)° (right) branching product. 

Similarly, for the study of Ar(3872)-like particles, we 
performed a mass scan inside the fitting region, i.e., a 
sequence of fits similar to the ones described above but 
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with a mass floating in 20 MeV/c^-wide mass bins. No 
significant signal is found in any of the studied bins. 

A similar analysis was performed by BaBar , where 
the upper limit on the product cr (77 —AT (3872)) x 
S(A'(3872) —>• r]cTr~^Tr~) was set. 

In conclusion, we report a study of the following B 
decays to final states with rjc- —>■ 

K^rjcOJ, B^ —>■ K^tJcT] and B^ —>■ K^rjcir^. 

We first study these B decays without intermedi¬ 
ate resonances and set 90% C.L. upper limits on their 
branching fractions: B{B^ —>■ K^'qcTT'^TT~) < 3.9 x 10“^, 
B{B^ < 5.3 X 10-^ B{B^ ->■ K^r]^r]) < 

2.2 X 10-^ and B{B^ < 6.2 x 10"^ 

We then assume that the decays proceed through in¬ 
termediate molecular states similar to the exotic X (3872) 
particle, such as (3872) = D°D*°-p°D*°,_X{3730) = 
-b D°D° and A:(4014) = D*°D*^ + D*°D*° and 
search for these states, setting 90% C.L. upper limits 
on the following product branching fractions: B{B^ —>■ 
K^X) X B{X -)■ r]cTr+Tr-) < 3.0 x and B{B^ 
K^X) X B{X ricLo) < 6.9 x lO’^ for X = Xi(3872), 
B{B^ K^X) X B{X rjcTj) < 4.6 x 10-^ for 

X = X(3730) and < 3.9 x for X = X(4014), 
B{B^ X±X) X B{X ? 7 c 7 rO) ^ 5 7 ^ 10“® for 

X = X(3730) and < 1.2 x 10"® for X = X(4014). 

We search for the neutral partners of the recently ob¬ 
served exotic states Z(3900)^ and Z(4020)^ and set 90% 
C.L. upper limits on the product branching fractions 
B{B^ K^Z) X B{Z 77c7r+7r-) of 4.7 x for 
Z = Z(3900)° and 1.6 x 10"^ for Z = Z(4020)°. 

We set 90% C.L. upper limits on the following B de¬ 
cays involving the state X(3915), whose origin is still un¬ 
known: B{B^ X±X(3915)) X B(X(3915) ^ 7 ^??) < 

3.3 X 10-5 and B{B^ X±X(3915)) x S(X(3915) ^ 
7?c7r°) < 1.8 X 10-5. 

There are no theoretical predictions for the decay 
branching fractions of D^D*^ molecular states similar 
to X(3872). In this paper, we obtain an upper limit on 
the product branching fraction ,6(5=*= —X=*=Xi(3872)) x 
6(Xi(3872) —>• ?7c7r+7r-), which is of the same order as 
the product branching fraction 6(5=*= —>■ X=*=X(3872)) x 
6(X(3872) —>• J/ijjTY^Tr~) measured in Ref. [s^. A sim¬ 
ilar situation is observed with the upper limit on the 
6(5=t ^ X=tX(3915)) X 6(X(3915) ^ r]^uj) and the 
value of 6(5=5 ^ iy±x(3915)) x 6(X(3915) J/'tpuj) 
obtained in Ref. [l^. A more copious data set expected 
from the upcoming Belle II experiment can pro¬ 
vide an opportunity to determine the ratios of the decay 
branching fractions of the exotic states described above. 
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